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ABSTRACT 

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data obtained as part of the Science Demonstration Phase (SDP). The 
data cover a central 4x4 sq deg region of the cluster We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 /im 
luminosityfFunctions (LPs) for optically bright galaxies that are selected at 500 //m and detected in all bands. We compare these LPs with those 
previously derived using IRAS, BLAST and Herschel-ATLAS data. The Virgo Cluster LPs do not have the large numbers of faint galaxies or 
examples of very luminous galaxies seen previously in surveys covering less dense environments. 
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1. Introduction 

The Virgo Cluster provides us with a unique opportunity to study 
in detail a large number of galaxies in the cluster environment. 
Virgo is probably the most studied cluster of galaxies because 
of its proximity to the Milky Way - it lies at a distance of ~17 
Mpc (Gavazzi et al., 1999 and references therein) with a mean 
velocity of -1094 km s"' (BinggeH et al., 1987). It is an Abell 
richness Class I cluster containing ~2000 optically catalogued 
galaxies (Virgo Cluster Catalogue (VCC), Binggeli et al., 1985). 
The Virgo Cluster virial mass is » 2.5 x 10'"^ M© (Girardi et al., 
1998; Rines & Diaferio 2006) as expected for a cluster as rich 
as this. The cluster contains large amounts of X-ray emitting gas 
(Boehringer et al., 1994) and shows clear evidence for both sub- 
structure and non-virialised motions (Gavazzi et al., 1999). From 
previous observations of galaxy clusters we know that this envi- 
ronment can have a significant effect on the properties of galax- 
ies. The morphological mix is quite different compared to the 
field and many galaxies may be deficient in atomic gas (Haynes 
& Giovanelli, 1984). This deficiency is probably responsible for 
the lower star- formation rates (Lewis et al., 2002, Gomez et al., 
2003), the truncation of star-forming discs (Boselli & Gavazzi, 
2006) and the observed higher metal content of cluster galaxies 
(Skillman et al., 1996). Being so nearby, the Virgo Cluster has 
previously been studied in detail at many wavelengths and so 
there is a wealth of data available, to which we are now able to 
add the FIR. 

This is the first of a series of seven papers, based on HeViCS 
SDP data, in which we discuss: how the cluster environment 
truncates the dust discs of spiral galaxies (paper II, Cortese et 
al., 2010), the dust life-time in early-type galaxies (paper III, 
Clemens et al., 2010), the spiral galaxy dust surface density and 
temperature distribution (paper IV, Smith et al., 2010), the prop- 
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erties of metal -poor star- forming dwarf galaxies (paper V, Grossi 
et al., 2010), the lack of thermal emission from the elliptical 
galaxy M87 (paper VI, Baes et al., 2010) and the far-infrared 
(FIR) detection of dwarf elliptical galaxies (paper VII, De Looze 
et al., 2010). A further paper discusses the spectral energy dis- 
tributions of HeViCS galaxies together with others observed as 
part of the Herschel Reference Survey (Boselli et al., 2010). In 
this paper we briefly discuss the HeViCS project, describe the 
data used in the above papers and compare the Virgo FIR LFs 
with those previously derived. 

2. The HeViCS 

The HeViCS is an approved Open Time Key Project for the 
ESA Herschel Space Observatory (Pilbratt et al., 2010). It has 
been awarded 286 hours of observing time in parallel mode with 
PACS (Poglitsch et al., 2010) at 100 and 160 fim and SPIRE 
(Griffin et al., 2010) at 250, 350 and 500 yum. We will map four 
4x4 sq deg regions of the cluster (Fig. 1) down to the 250 yum 
confusion limit of a; 1 MJy sr"' . From within these four regions 
we expect to detect about 200 galaxies in all bands. At the dis- 
tance of the Virgo Cluster a typical 50 kpc diameter galaxy sub- 
tends an angle of about 10'. Herschel with resolution of order 
10"is thus able to resolve structure and we can study the detailed 
FIR properties of these galaxies for the first time. 

The primary HeViCS science goals include the detection 
of dust in the inter-galactic medium, the extent of cold dust 
in the outskirts of galaxies, the FIR LFs, the complete spec- 
tral energy distributions (SEDs) of galaxies, the dust content 
of dwarf ellipticals and irregulars and a detailed analysis of 
the dust content of early type galaxies (for further details see 
http://www.hevics.org) . 

3. The data 

We have obtained SDP data for one field at the centre of the 
Virgo Cluster (Fig. 2) using the SPIRE/PACS parallel scan- 
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Fig. 1. Virgo Cluster. The solid black line indicates the area cov- 
ered by the VCC, with each VCC galaxy marked by a small 
black circle. The ROSAT X-ray contours are marked in blue. 
The area to be observed by HeViCS is indicated by the four red 
boxes. The SDP data described here are from the field centred 
on M87. 



map mode with nominal detector settings and fast scan speed 
(60"/sec). Two cross-linked scans were carried out, which is just 
one quarter of the exposure time of the completed survey. 

PACS data reduction was carried with the standard pipeline 
for both the 100 and 160 jjm channels. Dead and saturated pix- 
els were masked. Deglitching was performed in two steps, us- 
ing the standard multi-wavelength median transform deglitcher, 
and one based on sigma-clipping. Bright sources were masked 
before a high-pass filter was used to reduce 1// noise. Finally, 
the two orthogonal scans were combined and maps made using 
the naive map-maker. The full width half maximum (FWHM) 
beam sizes are 7.0"xl2.7"and 11.6"xl5.7"with pixel sizes of 
3.2"and 6.4"for the 100 and 160 //m channels respectively. 

For SPIRE the data were reduced up to Level 1 as described 
in Pohlen et al. (2010) except for the deglitching. Instead of the 
default wavelet deglitching module we used the newer kappa- 
sigma de-glitcher. Rather than the default median baseline sub- 
traction we used a simple robust linear fit with outlier rejection to 
all data points in the timeline. This gave a superior result in terms 
of residual striping. Both scans were then combined to make our 
final maps using the naive mapper provided. The FWHM of the 
SPIRE beams are 18.1", 25.2", and 36.9"with pixel sizes of 6", 
10", and 14"at 250, 350, and 500 yum, respectively. The uncer- 
tainty in the flux calibration is of the order of 15% (Swinyard et 
al., 2010). 

The final data used by the seven He ViCS papers have a 1 cr 
noise determined from the whole of each image of ~8.6, 4.9, 
1.3, 0.6 and 0.3 MJy sr"' at 100, 160, 250, 350 and 500 jum 
respectively. 



4. The luminosity functions 

To obtain LFs we chose to carry out the initial object selection 
at 500 yum because it is the least explored part of the spectrum. 




Fig. 2. 4" X 4" SPIRE 250 yum SDP data. Circles mark the po- 
sitions of the galaxies selected with 500 fim fluxes > 0.2 Jy and 
an 1.50-500 angular diameter > 1.4'. The cluster centre is marked 
by the position of M87. Most prominent is the extent of galactic 
ciiTus even at this high Galactic latitude (b « 75°). 



it has the lowest resolution, and most galaxies will produce their 
least flux in this band, guaranteeing a detection in all five bands. 
For the LF analysis the data were smoothed and re-gridded to 
the 500 yum resolution and pixel scale. The final Icr noise in the 
smoothed and re-gridded images are ~2.2, 2.9, 1.2, 0.7 and 0.3 
MJy sr-i at 100, 160, 250, 350 and 500 fim respectively. 

To produce an objectively selected sample we used the auto- 
matic image detection algorithm SExtractor (Bertin & Arnouts, 
1996). We detected at 500 fj.m and then used the same apertures 
at the other wavelengths. To minimise background contamina- 
tion each object had to have more than 30 connected pixels at 
1 .5cr5oo or above. The result is a 500 yum flux limit of ~0.2 Jy for 
sources with diameter larger than 1.4' and 97 detected objects. 
Each object was then checked for correspondence with a known 
Virgo Cluster galaxy. The final list consists of 24 Virgo Cluster 
objects, just 10% of the 247 VCC galaxies in the field that are 
listed by GOLDMINE as definite cluster members (Gavazzi et 
al., 2003). By far the majority of discarded objects appear to be 
associated with extended galactic cirrus emission. The exception 
are nine detections that have images that appear to be galaxies 
but are not in the VCC. In each case they are close to (~ 10") a 
faint SDSS galaxy with no redshift, so provisionally we assume 
these to be objects in the background. 

We also visually inspected the positions of all confirmed 
VCC galaxies in the 500 yum image to see if any other galax- 
ies could be detected. One hundred and fifty-eight (64%) of the 
VCC galaxies in this field are classified dE, one of which is a 
low signal-to-noise detection at 500 yum (De Looze et al., 2010). 
Twenty-eight (11%) are Sm/Im/BCD/dliT four of which are de- 
tected at 500 yum, but below our threshold (Grossi et al., 2010). 
Thirty-seven (15%) are other spiral and early type galaxies that 



Davies et al.,: The Herschel Virgo Cluster Survey: I. Luminosity functions 



with the exception of three are again not detected at 500 j-im. 
These three (1%) are clearly detected, but fail by their isophotal 
size. The positions within the cluster of the twenty-four 500 fim 
selected galaxies is shown in Fig. 2. 

The LFs at each wavelength are shown in Fig. 3. These FIR 
LFs are quite different compared to the optical, where large num- 
bers of dwarf galaxies are found at the faint-end. They are sim- 
ilar to the HI mass function, where there is a suggestion of a 
turnover at low masses (Davies et al., 2004). The comparatively 
low numbers at the faint end cannot be due to incompleteness 
because for this bright galaxy sample even the two faintest FIR 
sources have a signal-to-noise at 500 jum of ~20. It is clear that 
there are many low-luminosity galaxies in the cluster that are de- 
tected in the optical, but their lack of gas and dust makes them 
difficult to detect at 21cm and in the FIR. 

For comparison we include some previously determined LFs 
(Fig. 3). Rowan-Robinson et al. (1987) derived an IRAS 100 yum 
LF for optically identified galaxies (over a wide range of en- 
vironments) with photographic magnitudes less than 14.5 and 
an available redshift (blue dashed line Fig. 3a). This is roughly 
comparable with our Virgo sample in the sense that the galax- 
ies all have to have an optical identification and a distance. The 
Rowan-Robinson et al., LF cuts off quite sharply at the bright 

but then through to about 4 x 10^' 



end at about lO^^ W Hz"' sr"' 



W Hz sr it follows a power law slope of a;-l. Although the 
LF normalisation is arbitrary it is clear that the Rowan-Robinson 
et al., LF predicts ~ 10^ more galaxies at the faint end than at 
the bright end - even with the small numbers this is inconsistent 
with the Virgo data. We suggest that this is most likely due to 
processes in the cluster that remove gas and dust, i.e. tidal and 
ram pressure stripping (Doyon & Joseph, 1989). These processes 
are more effective for lower mass objects. The Rowan-Robinson 
et al., data also include galaxies that are one order of magnitude 
brighter than the brightest galaxies found in the cluster. Although 
in this sample there are very low numbers it has previously been 
noted that clusters in general do not have the very bright FIR 
sources seen in the field (Bicay & Giovanelli, 1987). Also in- 
cluded in Fig. 3a (light blue stars) are the IRAS 100 //m data, i.e. 
the LF for all VCC galaxies in the field with an IRAS 100 yum 
detection (31 galaxies). Although subject to very different selec- 
tion criteria the two data sets are broadly consistent with each 
other. 

In Fig. 3b we show the derived LF at 160 fim. We have pre- 
viously chosen to select our sample in a way that guarantees a 
detection in each of the five bands, but we are aware that obser- 
vational selection may affect any inferences made. To investigate 
this we used the programme SExtractor to detect objects down 
to adiameterof 30"at a 1.5cri5() isophote on the unsmoothed 160 
fj.m data. Note that much of the Galactic ciiTus is removed in the 
RA^CS data by the standard pipeline 1// noise removal. The re- 
sult is 1 12 detections of which 26 correspond with known VCC 
galaxies and are thus confirmed cluster members. These 26 are 
plotted as blue stars on Fig 3b, they are consistent with our pre- 
vious sample. The complete (112) sample (assuming they all lie 
at the distance of Virgo) are also plotted as green triangles, the 
dashed line indicating a signal-to-noise of 5. The numbers rise 
steeply, but with a slope less than expected for the background 
counts (slope of -1.5). The nature of these galaxies and similar 
ones detected in other bands will be investigated in subsequent 
papers. 

Bales et al., (2009) have used data from BLAST to derive 
250, 350 and 500 fim LFs for galaxies at different redshifts. On 
Figs 3 c, d, e we compare the lowest redshift (0 < z < 0.2) 
LFs (blue dashed line) with the Virgo data (arbitray normali- 
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Fig. 3. Luminosity functions at 100 to 500 fim for the Virgo 
Cluster bright galaxy sample (red crosses) with V(^) counting 
errors. The various other lines and points are described in the 
text. 
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sation). There are some very luminous galaxies (for example 
L250 > 10^"* W Hz"' sr"') in the BLAST data that are not 
found in Virgo, though their space densities are low at ~ 10"^ 
Mpc"^. The BLAST LFs do not go beyond the peaks in the Virgo 
Cluster LFs and so it is difficult to say what the predicted num- 
bers are at lower luminosities, but Bales et al. also show that the 
BLAST LFs are consistent with an extrapolation of the IRAS 
and SCUBA LFs, which do have a rising power law slope at the 
faint end. So our expectation is that we will again find that there 
is a lack of faint sources in Virgo when compared to galaxies 
in the general field. The solid blue line on Fig 3c is the recent 
Herschel-ATLAS (Bales et al., 2010) determination of the local 
250 yum LF (Dye et al., 2010) which is roughly consistent with 
the BLAST data. We look forward to a future comparison with 
other galaxy clusters at these wavelengths. 

Herschel large area surveys like ATLAS will eventually de- 
tect many more galaxy clusters. For example in the ATLAS 
fields there is the prominent and well studied Coma Cluster 
and ~20 other Abell clusters with z<0. 1 . The importance of the 
above Virgo data is that we can be quite sure that all of the galax- 
ies are cluster members - this is not so for surveys of more distant 
clusters. Bven so, with about 20 ATLAS clusters one would ex- 
pect that the cluster LF will eventually be well measured, even 
though the greater distances may inhibit this analysis. At the 
distance of the Coma Cluster, one of the nearest clusters in the 
ATLAS fields, (95 Mpc for Ho =73 km s"' Mpc"') our faintest 
500 fivn source (0.2 Jy at 17 Mpc) would have a flux of ~6 mJy. 
This is below the ATLAS 500 //m flux limit of 53 mJy. The 
brightest source detectable by ATLAS at the distance of Coma 
would have a luminosity of 5 x 10^' W Hz"' sr"', as indicated 
by the dashed line on Fig. 3(e). Thus almost 80% of the Virgo 
Cluster galaxies in this sample would not be detected at the dis- 
tance of Coma. The brightest source in Virgo with a 500 fim 
flux of ~12 Jy would fail to be detected beyond a distance of 
255 Mpc or z=0.06. Most of the clusters in the ATLAS fields 
lie within z=0.1, but at or beyond the distance of Coma, so our 
Virgo Cluster data predicts that only the very brightest cluster 
galaxies will be detected in the ATLAS. 

Using a single temperature modified blackbody (Ji = 2.0) 
and an emissivity of 0. 19 m^ kg"' at 350 fim (Draine, 2003) we 
have fitted the SBDs of the 22 late-type galaxies in the sample. 
Dust masses and temperatures within these 500 jum determined 
apertures are in the range of 10^^"^^ M© and 15-24K respec- 
tively. This illustrates the potential sensitivity of the full depth 
survey to low dust masses (< 10^ Mq) and the existence of a 
significant cold dust component (T< 20K). 

5. Conclusions 

We have described the HeViCS and our first science data from 
PACS and SPIRB. The data were used to produce 100, 160, 250, 
350 and 500 fim LFs for a sample of optically bright galaxies in 
the central region of the Virgo Cluster. Unlike the global optical 
LF and those previously derived in the FIR these functions do 
not appear to have a steep power law dependence at the faint 
end. There are relatively few faint FIR sources that can be 
associated with confirmed Virgo Cluster members compared to 
what is found for field galaxies and there are no examples of 
very luminous FIR sources. 
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